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CALCULUS OF VARIATIONS 



CHAPTER I 

INTRODUCTION 

1. The Calculus of Variations owed its origin to the attempt to 
solve a very interesting and rather narrow class of problems in 
Maxima and Minima, in which it is required to find the form of a 
function such that the definite integral of an expression involving 
that function and its derivative shall be a maximum or a mini- 
mum. 

Let us consider three simple examples: The Shortest Line, The 
Curve of Quickest Descent, and The Minimum Surface of Revolu- 
tion. 

(a) The Shortest Line. Let it be required to find the equation of 
the shortest plane curve joining two given points. 

We shall use rectangular coordinates in the plane in question 
taking one of the points as the origin. Call the coordinates of the 
second point Xi_, yi. 

If !/ = fix) is a curve through (0, 0) and (xi, yi) and I is the 
length of the arc between the points, obviously 

I = /"v dx^ + dy^ 
a 

or 7 = fv 1 + y'^dx, (1) 

and we w^sh to determine the form of the function / so that this 
integral shall be a minimum. 

(b) The Curve of Quickest Descent. Let it be required to find the 
form of a smooth curve lying in a vertical plane and joining two 
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2 CALCULUS OF VARIATIONS 

given points, down which a particle starting from rest will sHde 
under gravity from the first point to the second in the least possible 
time. 

We shall use rectangular axes in the vertical plane taking the 
higher point as the origin and taking the axis of X downward. Call 
the coordinates of the second point a^i, j/i. 




Let y — f{x) be a curve through (0, 0) and (xi, j/i) and use the 

ds 
well-known fact that tt' the velocity of the moving particle at any 

time, is ■^2gx. 



yl2(,x 



J ■<2n: 



/V 1 + y''^dx, 



and the form of the function / is to be determined so that this 
integral shall be a minimum. 

(c) The Minimum Surface of Revolution. Given two points and 
a line which are co-planar, let it be required to find the form of a 
curve terminated by the two points and lying in the plane which, 
by its revolution about the given line, shall generate a surface of 
the least possible area. 
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INTRODUCTION 3 

Take the line as the axis of X and use an axis of Y through one of 
the points. Call the coordinates of the points 0, ye, and Xi, y,. Let 
y = /(^) be a curve through (0, yo) and (3:1, j/i). 

If iS is the area of the surface of revolution generated by the 
curve, 

S = 2ir fyds = 2ir fy V I + y'Mx. 



-f^ 



and we wish to determine the form of the function / so that I shall 
be a minimum. 

2. The three problems just considered are special cases of what 
we shall call our fundamental problem which is, to determine the 

form of the function /so that if J/ = fix), j ^ (x, y,y') dxshaW be a 

maximum or a minimum; if being a given function and Xo and Xt 
being given constants, as are j/o and j/i, the corresponding values 
ofy. 

In ordinary problems in maxima and minima y = f{x) is a given 
function and we wish to find a value, xa, of x for which y is greater, if 
we seek a maximum, less, if we seek a minimum, than for neighbor- 
ing values of x; that is, for values of x differing from Xo by a suffi- 
ciently small amount whether that amount is positive or negative. 

In our new problems, to speak in geometrical language, we have 
to find the form of a curve for which our integral, I, is greater or 
less than for any neighboring curve having the same end-points. 

3. Let us now attack our first problem, that of the shortest line. 
We have to find the form of the function / so that if 



'f 



j ^l + y'-'dx 
I shall be a minimum when y = f{x). v. Art. 1 (a). 
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4 CALCULUS OF VARIATIONS 

Let y = F{x) be any other continuous curve joining the given 
points, and let ri{x) = F(x) — f{x). 

Then y = f(x) + -ni^) is our curve y = F(x). Consider the 
curve y = f{x) + a?j(x) where a is a parameter independent of x. 

Y 




For any particular value of a the curve y = f(x) + aijix) is one of 
a family of curves including y = f{x) f or o = and y = F{x) for 
a = 1. 

By talcing a sufficiently small value for o we can make aij(x) less 
in absolute value for that and all leas values of a, and for all values 
of X between and Xi, than any previously chosen quantity ^; and 
for such values of a the curve y = fix) + aijix) is said to be a 
curve in the neighborhood of y = f{x). 

liy = /(*) andi/ = F(x) are given, 7(n), the 7 for any one of our 

curves y = f(x) + aijix), is f Vl + (i/' + ar]'(x)ydx, and 1(a) is a 
function of a only. o 

A necessary condition that 1(a) should be a minimum when 

a = is well known to be tliat -r- 7(a) should be zero when a = 0, 

This condition we shall express as /'(O) = 0. 
Since the limits and Xi are constants 

/'(«) -/^ ^T+ly^+^'^^dx 



and /'(O) = / ,_^ yi'ix)<tx. 

J V 1 + y'2 
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INTRODUCTION 5 

Integrating by parts 

L V 1 + 3/'2 J J^"^ Vl + y^ 

J dx Vl+y'^ 

since )j(x) vanishes when « = and when x = Xi. 

A necessary condition that 1(0) shall be less than /(a) for some 
value of o and for all less values of a no matter what F{x) may be, 
in which case the length of the curve y = f{x) is less than that of 
any neighboring curve, is that /'(O) = independently of 7i{x), 

i. e. that / -r — ■ - ri(x)dx = no matter what the form of the 
J da; V 1 + y'^ 

arbitrary function 7i(x). 

This condition will be satisfied if and only if -r- , = ; 

and we thus are led to a differential equation of the second order 
between y and x. 

Its solution will express y as a function of x involving two arbi- 
trary constants. , 



whence y' = K, a constant; 

y = Kx + L. 
The required curve is to pass through (0, 0) and {xi, yi) and thus 
we are able to determine K and L. 

L = Q, K=^- 

X\ 

Hence y ^ ^ x; and our curve is a straight hne through the given 
points. ' 

If certain other conditions, depending on the fact that when 1(a) 

is a minimmn /"(o) must be positive, are satisfied, y = —x must 
be the required shortest line. ^ 

As our necessary conditions gave us but a single solution it is 
clear that if there is any shortest line it must be our line y = — x. 



y Google 



6 CALCULUS OF VARIATIONS 

We may note in passing that in simplifying J'(«») by integration 
by parts we tacitly assumed that j{x) and F{x) were continuous 
and had continuous first derivatives over the range of integration. 

4. We can deal with the general problem formulated in Art. 2 
precisely as we have dealt with the shortest line problem. 

Let it be required to determine the form of the function /so that 

y = f(x) shall make j * (x, y, y') dx a maximum or a minimum; 

given that y = yo when x = Xa and y = yi when x = Xi. 




As in the last article let tj(x) = F(x) — fix), and consider the 
family of curves 

y = fix) + ariix). 

Let 1(a) = / 4>(x, y + aij{x), y' + aji'ix))dx'. 

I'{a) = / ^ 4>(x, y + av(x), y' + aTi'{x))dx. 

Let u ^ y + aij(x), v = y' + av'ix). 

du , , dv ,, , 

/'(a) = / T- fPix, U, v) dx 
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INTRODUCTION 



and /'(O) must be made equal to zero. 
Integrating by parts, 

7'(0)=/|,M.. + [*,wf-/||,W*, 

since v{^) vanishes when a: = a^ and when x = xi. 

/'(O) must be zero independently of the form of ij((r) therefore 



dx i 



-- 0, (2) 



and as before we are led to a differential equation of the second 
order which if solved gives 3/ as a function of x involving two arbi- 
trary constaits which must be determined from the facts that 
y = ya when a; = a^ and y = yi when x = Xi. This differential 
equation is known as Lagrange's Equation and it is a necessary con- 
dition that / should be a maximum or a minimum. 

Any particular solution of Lagrange's Equation is ctilled an 
extremal, and if the given problem has a solution it is that extremal 
which passes through the given end-points., 

If cj is a function of y' only Lagrange's Equation becomes 

I -*«-«■ 

Whence y' = k, 

a constant, and the extremals are straight lines; and therefore the 
required solution is the straight line through the given end-points 
as in Art. 3. 

The problems of Art. 1 (6) and (c) can be solved by substituting 
in Lagrange's Equation the appropriate value of ^ and then solving 
the resulting equation. 
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\ CALCULUS OF VARIATIONS - 

For the curve of quickest descent 

= ~- -ll + y" V. Art. 1 (b) (2). 

d4> 1 if 



h' ^fx V 1 + j/'3 
's Equation becomes 

f -■= ^=«=- 0. (3) 

die 4x V 1 + y'2 

For (/te minimum surface 

i, = V^TT7' I-. Art. 1 (c) (3). 

8^ ^ yy' 
<ni' V 1 + !/"' 

Lagrange's Equation becomes 

-. , I .... i ra' _ o /<> 



We shall reserve the solving of (3) and (4) for a later article. 
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CHAPTER, II 

VARIATIONS. NOTATION AND NOMENCLATURE. 
ILLUSTRATIVE PROBLEMS 

5. If y = f{x) and y = F{x) are the equations of two curves and 
ij{a;) = F{x) — f{x), ii(x) can be regarded as the increment pro- 
duced in y by changing f{x) into F{x) without changing x. We 
shall call it the variation of y and represent it by Sy. It is a function 
of X and usually a wholly arbitrary function of x. 




n'(x) = -J- Sy and is easily shown to be the increment produced 

in y' by changing /(x) into F{x), x being held fast, since when y 
becomes y + t){x) y' becomes y' + ri'(x) and is changed by ij'Ca;}. 
We shall call this increment the variation of y' and represent it by 
5y', Like 5u it is a function of x. As we have seen 



;«». 



(1) 



If y is increased by Sy^ 



is an approximation to the increment produced in ^(j/), the close- 
ness of the approximation depending on the magnitude of &y. If 
Sy is infinitesimal this approximation differs from the actual incre- 
i^icr.* by an infinitesimal of higher order than dy. 
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CALCULUS OF VARIATIONS 

II call this approximate increment the variation of ^(t/) and 
represent it by Sipiy) ; so that 



rr an approximation to the increment produced in 
(l>(y, y') by giving y the increment &y and y' the increment &y'; the 
closeness of the approximation depending on the magnitudes of &y 
and hy'. We shall call it the variation of i^(y, y') and represent it by 
b^iy, y'), so that 

As ail our variations are supposed to be caused by changing the 
form of the function f my = J{x) without changing x, 



If in (2) and (3) we replace the symbol 3 by the symbol d we 
obtain the familiar formulas for d4>. Hence variations are calcu- 
lated by the use of the familiar formulas and processes of the 
Differential Calculus. Indeed if j/ is a function, ff of x, d4> and Scj) 
are approximate increments of precisely the same type but dif- 
ferently caused, rf*, by changing x without changing the form of 
the function and S<i>, by changing the form of the function without 
changing x. 

With our new notation our necessary condition that /, the 

<j>{x, y, y') dx, shall be a maximum or a minimum can be written 
60(a:, y, y')dx = 0, (5) 



/• 



V. Art. 4 (1). 

If we define the variation of a definite integral as the integral of 
the variation of the integrand, so that 
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ILLUSTRATIVE PROBLEMS 11 

S j 4>.dx = I 5il>.dx, our condition (5) 

may be abridged into 

57 = 0. (6) 

6. To illustrate the ordinary technique of the Calculus of Varia- 
tions we shall now work out at length the problem of Art. 1 (fe) ; 
the bracfiistochrone or curve of quickest descent. 

We wish to make / a minimum where (v. Art. 1 {(»)) 



I 

Wo must make Si = 0. 



^ \S V 1 + w'^ «^ ■ 

^ 

Integrating by parts, 

LVS V 1 + y'2 Jy J^<i^l^x VI + y'^ 

Since our end points are fixed Sy — when a: = and when 
X = a:i. Therefore ^^ 

5l=~ff\^-^]8ydx. 
J dxLvi V 1 + y'U 

As 3/ must be equal to zero no matter what the value of our arbi- 
trary function Sij we must have 

dxl-j^ V 1 + y'sj 
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12 CALCULUS OF VARIATIONS 

and (1) is our differential equation for y, i. e. tlie differential 
equation of our required curve, v. Art. 4 (3), 
(1) is easily integrated. We have 




y/l. 





y 


V 


xit 




- - V 2S 
when X 


= 0, 

+ a 


+ a vers- 

c . 

ver.,-.i, 





+ (" 


Since 


ax - 
-0 

V - 






and 


V2<ii - X' 





(2) 

where a is to be determined so that the curve shall pass through the 
point {xi, yi), is our required curve of quickest descent. 

(2) is the equation of a cycloid having its base in the axis of Y, 
its cusp at the origin, and lying on the positive side of the axis of Y. 




Our solution, then, is an inverted cycloid with a horizontal base and 
it has a cusp at the higher of the given points. 

Our treatment has assumed that in the required curve y is a 
single-valued function of x. Our solution is, therefore, apparently 
valid only if a cycloid suitably situated can be drawn so that 
(xi, yi) shall lie between the origin and the vertex of the curve. As 
a matter of fact this limitation is unnecessary as may be shown by 
interchanging the axes of X and Y and working the problem anew. 
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ILLUSTRATIVE PROBLEMS 



If this is dono / = / ^r-^ dx and it is not difficult to make 



r v 1 + y'-- 
•I v^ 



SI = 0. The resulting differential equation 
2^" + 1 + I/'* = 
can be solved by the method employed in the following section and 
we find as the equation of our required curve 

X = - ■J2ay - y^ + avers"!-- (3) 

Throughout the whole arch of the cycloid (3) y is a single valued 
function of x, and {xi, j/O is, therefore, not restricted to lying in the 
first half of the arch. 

In any concrete problem the numerical determination of a is a 
matter of some difficulty. This, however, may be avoided by 
Newton'a famous solution for the curve of quickest descent from 
AtoB. 

' ' Through A on the horizontal line A C as a base draw any cycloid 
cutting the straight line AB in Z). Through A on the line AC draw 



a second cycloid whose base and altitude are to the base and alti- 
tude ofthe first as AB is to AD; it will pass through S and will be 
the curve required." 

7. As another example we shall take the problem of Art. 1 (c), 
the minimum surf ace of revolution. Here {v. Art. 1 (c), (3)), 
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•/ L V 1 + y'^ J 



Integrating by parts, 

<^ ^1 + y 



f/[,T-— ^-JS^ 



-f[^^--i^h^- 



To make 5/ = we must make 
d 



-- 0, (1) 



and this is the differential equation of our required curve, v. 
Art. 4 (4). 

Performing the differentiation indicated in (1) and reducing, the 
equation becomes 

yy" ~1- y'^ = 0. (2) 

Equation (2) can be integrated by a regulation device as followa: 
Let s = y'. 

„ _dz _dzdy _ ,dz _ dz 





^ dx dy dx " dy l 


(2) becomes 


».| = i + .- 




rT7W-f- 




ilog(l+#) -logij-C. 




^1+-' „ i, 



where a is an arbitrary constant. 
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ILLUSTRATIVE PROBLEMS 



j/ = «coshr + M. (3) 

(3) is the equation of a catenary having the axis of a; as its direc- 
trix. The constants a and h must be determined so that (3) shall 
pa^ through the two given points (0, yo) and (xi, j/O- To calcu- 
late their iimnerical values in a concrete case is not easy but the 
curve required may be obtained mechanically by the following 



Insert smooth pegs, at the two given points, in the vertical plane 
through these points. Throw over the pegs a uniform flexible 
string so long that one end will hang below the given axis; and then 
pvdl down on the other end until the free end is drawn up to the 
axis. Then the portion of the string between the pegs will be the 
required catenary. 

For it is well known that if one end of a uniform flexible string is 
fastened and the other end is thrown over a smooth horizontal peg, 
the curved portion of the string when in equilibrium will be a 
catenary whose directrix is a horizontal hue through the hanging 
end of the string, 

8. Since any function of a single variable can be geometrically 
repr^ented by a graph by treating the independent variable as an 
abscissa and the corresponding value of the function as an ordi- 
nate, our theory covers the case where we are to make the definite 
integral of a given function of the independent variable, and of the 
dependent variable and its derivative a maximum or a minimum, 
no matter what the ordinary interpretation of the variables in 
question. 

For example let us work the shortest line problem (Art. 1 (a)) 
using polar coordinates. 
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16 CALCULUS OF VARIATIONS 

Let r = f(<p) be the required shortest plane curve connecting the 
■points (ro, <l»)) and (n, (jn). 



Since ds = V dj-^ + r 

and we must make 

51 = 0. 



SI = I aV^M^T^d^ 



P rSr + r'Sr' 






Integrating by parts, 

f , '' . irSr.d<l> = [ , '' 6r] - f^ - — '"-^ 6r# 
J Vr^ + )-'2 a<ii L V r^ + r'^ J ^ «* V j-^ + r'^ 

5r vanishes when (^ = ^o and when c^ = iji,. 

Therefore 5/ - f] , '" - -^ '' Isrrto. 

To make SI = we must make 

'' _ — ^ = 0- (1) 

and this is the diiierential equation of our required line. 

If we perform the indicated differentiation and reduce, (1) 
s the comparatively simple form 
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17 



We can clear of fractions and solve without difficulty, but if we 
happen to remember the Calculus formula for the curvature of a 
curve in polar coordinates 



[-(1)1 



we see that our differential equation (2) represonts a curve of zero 
curvature, i. e. a straight line. 

9. As another example let us find a geodesic line joining two 
given points on the surface of a sphere. 



The element of arc of a curve on the spherical surface is 
ds = V a?M^ + a* sin^ fld^^ 



and we may take J — f V 1 + sin^ fftp'^ c 
J V 1 + sin^ 50'^ 



J ^n~- 



Iiitegrating by parts and remembering that 3* vanishes at both 
limits 00 and Si, we have 



J ae VT+am^e^'^ 
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18 CALCULUS OF VARIATIONS 

Since 51 ^ whatever the form of &tj> 
sin^ dtp' 



Integi 


■ating, 
C 


VBC0 


<^evi + sin^e^'^ 


If we 


- C V 1 + sm= 9,^'^ 
C c)5c- e 


let z = 


■■e - C^ V 1 - C^ csc^ 
ctii and K^ = . ^^ 



0- (1) 



Vl -6'^- CcWt 



-I. 



*'" - - sin-'fc + a - - ain-' (KctiiS) + « 



where a is an arbitrary constant. 

sin {a — tp) = K ctn 0, 
sin a sin S cos ^ — cob a sin fl sin ^ = K cos 0. (2) 

If we multiply by a and change to rectangular coordinates we get 
ysina ~ z cos a — K,y, 

a plane through the centre of the sphere. The constants must be 
determined so that this plane shall pass through the given points 
and then the great circle in which it cuts the sphere is the required 
geodesic. 

Examples 

(1) Work the bracMstochrone problem taking the axis of Y verti- 
cally downward, v. Art, 6, page 12. 

(2) Integrate it" — 2r'^ — r^ — 0, the polar differential equation 
of the shortest line. v. Art. 8 (2). 

Ans. r cos (^ — n) = p; a straight line. 
"(3) Find the shortest line that can be drawn on the surface of a 
cyhnder of revolution and joining two given points on the 
surface; a geodesic. Suggestion. Use cylindrical coordinates 
a, tji, z, taking z as the independent variable. 
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Then / = fv 1 + aV¥s. 



Ans. The geodesic is the line of intersection of the cylin- 
der with the helicoid <> = fca + c; and is a helix. 
(4) Find a geodesic on a cone of revolution. 

Suggestion. Use spherical coordinates r, a, <p, taking r as 
the independent variable. 



/> 



Then / = / V 1 + r^ sin^ a-l>'Mr. 



Ans. The geodesic is the intersection of the cone with 

the surface r cos [(</) — (3) sin a] = c. 

(5)' A mountain is in the shape of a hemisphere and the velocity 

of a man walking upon it varies as his height above its base; 

show that his path of least time between two given points 

upon its surface hes in the vertical plane through the points. 

10. Isoperimelrical Problems. 

There is an important class of problems in which it is required 
to make I a maximum or a minimum while keeping constant the 
integral of a second given function of x, y, and y' taken between the 
same limits. 

(a) For example take the problem to find the closed plane curve 
of given perimeter and maximum area. 



We shall use polar coordinates and take the origin within the 
required curve. 

If >■ = /(*) is our curve, its area, A, and its perimeter, 5, are equal 

2ir 2i- 

respectively to / 1 i^d<j> and / V r^ + r'^d<p. 
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20 CALCULUS OF VARIATIONS 

We are to make A a maximum while S is kept equal to a given 
constant I. 

If X is a constant multiplier at present undetermined ^ + XS 
will be a maximum if our problem is correctly solved. 



Let, then, / = / [|-!^ + X Vr= + r"']d<j>, 

and proceed as usual. 



8I=f[r,r 



+ X- 




Performing the indicated differentiation and reducing we get 



as the differential equation of our required curve. 

The first member is the curvature of the curve, and as the cur- 
vature is ^, a constant, the required curve is a circle, v. Art. 8, 



As X is the radius of the circle its given perimeter I is SttX and 
our X proves to be s- ■ 
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ILLUSTRATIVE PROBLEMS 21 

(b) The ends of a uniform string are fastened at given points ; 
find the equation of the curve in which it must hang that its centre 
of gravity may be as low as possible. 

Let (xo, yo) and (a^i, yO be the end points, and let I be the length 
of the string. 



If y is the ordinate of the centre of gravity 

y = J j y V 1 + y'-dx. 

We must make i, 

y + \l, i. e. J J {y + \l) V 1 + y'Mx, 



a minimum. 



= U and let I = / (y + k)V 1 + y'^dx. 



Turning to Art. 7 we see that we have only to replace yhyy -}- k 
in the solution there given to get the solution of the present prob- 
lem, and that our required curve is y+K = a cosh f - + b), a 

catenary. Unlike the catenary in Art. 7 it does not generally have 
the axis of X for its directrix. 

In a numerical problem a and b would have to be determined so 
that the catenary would pass through (xo^.j/o) and (xi, j/i) and k 
would have to be found by the condition that 



rVT+Y"^<te = I. 



Problems like (a) and (6) are known as isoperimetrical problems 
and can all be handled by the device we have illustrated. 
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Examples 



(1) Find the form of a curve it the area bounded by an arc of 
given lei^th, the ordinates of its given end points, and the 
axis of X is to be a maximum, 

Ans. A circle. 

(2) Find the form of a curve if the sectorial area bounded by an 
arc of given length and the radii vectores of its end points is to 
be a maximum. 

Ans. A circle. 
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CHAPTER III 

PROBLEMS INVOLVIKG SEVERAL DEPENDENT 
VARMBLES 

11, It is easy to extend our theory to problems where instead of a 
single dependent variable y and its derivative there arc several varia^ 
bles to be determined as functions of the same independent varJa-> 

blc. That is, when I = f <i>i^, y,z,-- ■, y', z', ■ ■ •,)dx and 7 is to be 

made a maximum or a 

Let 






■ + ^V+^a.' + . 



and by the method of Art'. 4 it can be shown that we must make 
SI = 0. 



As an example let it be required to find the shortest c 
y plane joining (0, 0, 0) and (xi, y-,, Zi). 



* = ^1 dx^ + dy^ + dz^ = V 1 + j/'^ + z'^dx. 
' - fv 1 + y'= + z'Mx. 



^„ VI + y'^ + z'^ 



By integration by parts 



»--/t 



_dX^U + y'^ -I- 2'3 



<i": V 1 + !/'■ + 2 
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To make SI = we must make 

d v' . , d z' 

-T- ^ = By + :r t= &z 

vanish when Sy and &z are independent arbitrary functions of x. 
Hence 

|. .^^„_ . 

rfx V 1 4- i/'^ + 2'3 

and ^ — = 0. 

dx V 1 + j/'= + z'2 

These are a pair of simultaneous differential equations of the 
second order connecting y and s with x ; and can be solved without 

difficulty. We get t = d = 7* '^ straight line. 

12. Hamilton's Principle, called by Professor E. B. Wilson " the 
most fundamental and important single theorem in mathematical 
physics," and closely allied to the so-called principle ©f Least Action 
is that in the actual motion of any conservative system the time4ntegral 
of the sum of the Energy and the force function is a minimum. 

Formulating, 8 j [T + U]dt = 0, ii U is the force function and T 

the kinetic energy. 

We give a proof of Hamilton's Principle for a single moving 
particle. It can be extended without difficulty to a moving system. 

Let a particle of mass m moving under forces X, Y, Z, actually 
move from its initial to its final position in the time ti. Suppose 
that by the introduction of additional forces it had been made to 
move from its initial to its final position in the same time (i, but by a 
slightly different path. . Let (x, y, z) be its position at any time 
t in the actual motion and (a: + 5a;, j/ + Sy, 2 + hi) its position in 
the hypothetical motion after the lapse of the same time. If t is 
our independent variable we have from Mechanics mx" = X, 
my" = F, ms" = Z; and if the forces are conservative and V is 
the/orce /MJM^jon Xfis + F5j/ + Zhz = SU. 
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SE'V'ERAL DEPENDENT VARIABLES 
X5x = mx''Sx = m\-j-(x'Sx) — x'Sx'l = m-r,(x'Sx) 
5U = ^iiWSx + y'5y + z'bz] — 5T, 
where T ~ -^ [x'^ + y'^ + e'^] and is the Kinetic 



j S{T + U}dt = m[a;'5a! + y'fiy + z'Sz]_= 0, 



since the initial and iinal positions of the particle are the same in 
the actual motion and the hypothetical motion. 

As an illustration of the use of Hamilton's Principle we shall 
obtain the differential equations for planetary motion by its aid. 
If we use polar coordinates and take the sun as origin T, the energy 
of the planet, is -^ (r'^ + i^^'^). The force, F, exerted by the sun 
. is an attractive force directed toward the sun, proportional to the 
product of the masses and inversely proportional to the square of 
the distance. ^j^^ 

The work done in any motion of the planet is / Fdr 
and jFir=-fM«J%.'^ + C. 



By Hamilton's Principle 



m 



SV' 



•'&-j>' + r<i>"'Sr -—-Srldt^O, 
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Integrating by parts, 



Whence 
and 



Example 



A heavy cube containing a smooth spherical cavity rests on a 
smooth horizontal plane, and a heavy particle lies at the 
bottom of the cavity. The cube is given a horizontal velocity- 
Find the equations for the subsequent motion of the system. 
Suggestion. Let M be the mass of the cube and m the 
mass of the particle, and use as coordinates x, the distance 
the cube has moved, and 6, the angle a radius of the cavity 
drawn through the particle makes with the vertical. The 
horizontal, velocity of the particle is easily seen to be 
x' — a cos B6', its vertical velocity is d sin 8$'. 

T = i{M + m)*'^ + I'tc^^ - 2a cos Bx'8'. 
U = mga cos 6. 
The equations required are 



dt'' 



J \{M + 'm)x' — ma cos W\ = 

' - cos te') - sin 9 (x'd' - p) = 0, 

(M + m)x' ~ ma cos ee' = C 

oB" — cos to" + j7 sin = 0. 



13. If in the proof of Hamilton's Principle in the last article we 
do not replace XSx + YSy + ZSz by its equal dU we get 



/i 



[5T + X&x + YSy + 2 



and this equation is valid even if the forces are not conservative. It 
is to be noted that XSx + YSy + Z5z is the work that would be done 
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SEVERAL DEPENDENT VARIABLES 27 

by the actual forces, if the particle were displaced from its actual 
position at any time to it« corresponding position in the hypotheti- 
cal path. 

The differential equations in polar coordinates for the motion of 
a particle in a plane are easily obtained by the aid of our new equa- 
tion. Let the component of the force along the radius vector be R 
and perpendicular to the radius vector be *. If (r, i>) is the actual 
position of the particle at the time t the corresponding position in 
the hypothetical path ia (r + Sr, ij> + Sif). 

XSx + Y8y + ZSz 



f\^ 



j'lMrV 



I W^- nfr" + r4,' 



(d , 



(rVO - *J-)5^ di = 0. 



wr" — r-ti'^h^ 1 



(1) Find the differential equations for the motion of a spherical 
pendulum. 

Ans. 8" — sin 6 cos fl<^'^ + - sin = 



(2) Obtain the differential equations for the motion of a particle in 
space, using cylindrical coordinates. 
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(3) Obtain the differential equations for the motion of a particle ir 
space, using spherical coordinates. 

Ans. mlr" - r (fl'^ + sin^e.^'^)] = R 

~\m ^'^^'^ - '^ sin fl cos e^'^l = e 

m d . . , 
— ^-r -77 (r^ sin^ 60') = *. 
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CHAPTER IV 

MULTIPLE INTEGRALS 

14. If / is a multiple integral no serious additional difficulty is 
presented. We have still to make dl = 0. 

If I ~ j I <t>{x, y, z, p, qjdxdy, 

where P = ^=r 

'^ ox 

. & 

and q = ^-> 

dy 

and the integration is over a given area in the XY plane, let it be 
required to make / a maximum or a minimum by suitably determin- 
ir^ 3 as a function, /, of x and y, s being given for all values of x and 
y corresponding to points in the boundary of the area. 

Here &2 is the increment produced in z by charing the form of 
the function /, x and y being held fast, and is a function of x and y. 

If we let Bp = w- fe> and fig = ^fis then 5p and Sq are the incre- 
ments produced in p and q. 

2 = fix, y) + ahz is any one of a family of surfaces of which 
s = Six, y) is one. 

lia) = \j <l>ix, y,z + a5s, p + aSp, q + aSq)dxdy. 
= I f 54> dxdy = S j j ipdxdy = 37, 
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and, as before, the necessary condition that I shall bo a maximum 
or a minimum is Si = 0. 

15. Minimal Surface. As an example let us find the form of the 
surface of least area that can be bounded by a given closed curve 
not necessarily plane. 

I = f I Vl + p^ + (fdxdy, the integral being taken over the 

area bounded by the projection of the given curve on the XY plane. 



JJ ^l + 'p^ + q^ 



'If- 



" ^ ox cy . 

V 1 + p* + 9^ 



Integrating by parts, 



J V 1 + p* + 



= dydx. 



::dx = 



V 1 + p= + 9= 



J dx VT+p^ 



+ q' 



- dzdx. 



Since in our x-integration y is held fast our limits are the abscissas 
of the points on the projection of the given boundary which corre- 
spond to the value of y in question, and for these points 8s = 0. 
Therefore the term outside of the sign of integration vanishes and 



-5z 



J ^ll + p^ + q' J ^^ 4TTW 

it may be shown that 



J V 1 + pH 



J Si 



V 1 + p= + g^ 
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J J L^a: V 1 + p= + 2= ^y VI + p^ + f/J 

JJ (1 + p" + §')= 

3/ = when and only when 

[ \dy/ J 8x^ dx dy dxdy [ \dx} j g;/^ ^ ^^ ,^. 

and (I) is the differential equation of the required minimal surface. 
It is shown in works on Differential Geometry that the first mem- 
ber of (1) is the sum of the curvatures of any two mutually per- 
pendicular normal sections of the surface, and its half is called the 
m,ean curvature of the surface at any point. Therefore a minimal 
surface is a surface of zero mean curvature. 

Examples 

{!) Show that the mean curvature of a helicoid s = fc tan-i- is 
^ ' X 

(2) =^ cylindrical cup with a plane bottom of any shape and a given 
upper rim of any shape is to be capped so that the closed cup 
shall have a given volume and a minimum upper surface. 
Show that the cap must be a surface of constant mean curva- 
ture. 

16. The differential equation for small transverse oscillations of 
a stretched elastic string can now be obtained from Hamilton's 
Principle. 

Let I be the length and p the tension of the string in its position 
of equilibrium, and let the motion be so small that higher powers 
than the first of the slope of the string at any time may be neg- 
lected. Let us assume also that longitudinal motions and the 
change in tension due to the stretching of the string during its 
oscillation are negligible. 
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If m is the mass of a unit length of the string the kinetic energy 

T is equal to / -^ ( -^ ) dx. The force function U is the negative of 



the potential energy, V, of the string, and V is equal to the work 
that would be done by the tension, p, in restoring the string to its 
length when in equiUbrium: i.e. V = p (s — I), 





2 

SO that approximately V— / ^{a~) i^^- 



By Hamilton's Principle SI = 0, 
where I =J\t ~ V)dt. 

-//[i(ir-ii)']"- 









since &y = when t = and when ( = (. v. Art. 12. 
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dl - 

Hence 



MULTIPLE INTEGRALS 



'fib 






Obtain the differential equation for small transverse oscilla- 
tions of a stretched elastic membrane. 

««««*"■ 1' - '//v''i+(S" +(!)'"» -^ 

where A = j j dxdy is the area of the membrane at rest. 
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CHAPTER V 

VARIATION OF THE LIMITS. PRINCIPLE OF LEAST 
ACTION 

17. In our fundamental problem it is required to determine y as 

a function of x so that / i^(:t:, y, y')dx shall be a maximum or a 

minimum, xo, Xi, and the corresponding values ya and j/i being given, 
so that dy = when x = Xq and when x = Xi. 

If while keeping the range of integration fixed we remove the 
restriction that yo and j/i are given the old reasoning by which we 
established the necessary condition 

>■■« 

still holds good; but when we integrate by parts and get 
the last term no longer disappears, but becomes 



-/te 



where Syo and &y, are entirely arbitrary; so that dl = when and 
only when 



i] — [l^]-' 



and we must detennine the two arbitrary constants in the solution 
of the first equation by the aid of the second and third equations. 
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35 



The problem of sailing to windward is an interesting illustration. 
The speed of a sailing vessel is a function, more or less complicated, 
of the angle her course makes with the direction of the wind. It is 
required to find how she must sail to get a given distance, a, to 
windward in the least possible time. 

Let the starting point be taken as origin and an axis of X be 
drawn directly to windward. If j/ = /(x) is the required path y' is 




the tangent of the angle the course makes at any time with the 
wind's direction, and the velocity can be expressed in terms of y'. 




and I — I '■ — fT7~,T — dx. 





Since we must make 5l = 

_3_yT+Y 

dy' F{y') 



and 



ii'W-]-- 



(1) 

(2) 
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36 CALCULUS OF VARIATIONS 

(1) gives us y' = e, and the required path is a straight line. 

(p. Art.. 3.) 

(2) reduces to 

y'F(y') - (1 + y")F'(y') = 0, 

or since y' = c to cF(c) — (1 + c^)F'(c) = 0, 

which determines c, the constant slope of the required path. 

As a rough approximation of the speed law for a rather sluggish 
boat let 

v = F(,y') =/c(tan-'i/'~a). 




then, gives us the angle, 0, the course (known to sailors as full 
and by) should make with the direction of the wind; and as 6 
does not depend upon fc this course is independent of the strength 
of the wind. 

If the point aimed at lies within the angle formed by lines 
through the starting point making angles and — with the wind's 
direction it is reached by tacking, both tacks being steered full and 
by. If it lies outside of this angle it should be steered for directly. 

Examples 

(1) If a = I - 1 = 33° - , e = 62^°. 

(2) Find the curve of quickest descent from a given point to a given 
vertical hne. («, Art. 9, Ex. 1.) 

Ans. An inverted cycloid with a horizontal base, a cusp 
at the given point, and the vertex in the given Une. 

(3) Given two mutually perpendicular lines and a point in their 
plane; find the curve terminated by the point ajid one of the 
lines which by its revolution about the other line shall generate 
a surface of minimum area. (u. Art. 7.) 

Ans. A catenary having the axis of revolution as its 
directrix and having its vertex on the other given line. 



y Google 
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18. If 7 = j ij> (x', y, y') dx and we vary Xa and Xi as well as j/o 

and yi, so that the range of integration is no longer fixed, the prob- 
lem is more complicated. 

For the sake of simplicity we shall suppose that xo and yo are 
given but not Xi or j/i, and that we are to determine the curve 
y = f(x) so that / shall be less if we seek a minimum, greater if we 
seek a maximum than when y — J{x) is replaced by any neighboring 
curve and Xi is replaced by any value, xi + 5a:i, differing from Xi 
by a sufficiently small amount. 

With oxir ordinary notation y — f{x) + aSy is a neighboring 
curve if a is sufficiently small (v. Art. 3); and /(a), where 

/(a) == J-t>(x, y + aSy, y' + <xh/) dx, 

must be a maximum or a mmimum when a = 0. 
Let u = y + aSt/ and w = y' + a5y', 

ii + Qir, 

then I{a) = \4>{x,u,v)dx. 

and I'm =f(^^Sy + ^,&y'yx + [<t>ix,y,y')Ux, 

= j S<l,dx + [4>{x, y, y')]^Xi 
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To make 7'(0) = we have as usual Lagrange's Equation 



(1) 



dy dx dy' 
and we have the conditions 

y — Va when x = Xa, 

and [^, 6y] + [<l>{x, y, j/')]e^. = 0, (2) 

to determine the two arbitrary constants in the solution of 
Lagrange's Equation. 

If (Sy)i, the value Sy assumes when x = Xi, and 5;Ei are independ- 
ent the solution is generally impossible, but in many interesting 
problems (By)i depends upon 6xi, 

For instance suppose that the required curve y = f(x) which 
makes I a maximum or a minimum is to be terminated by a given 




point (xo, yo) and a given curve y = ^(a:), so that j/i = ipix^}. Let 
the second end point in the varied curve be {a^i + 5xi, yi ■+■ Syi). 

Then Syi = ^'(xi)5xi + e 

and (6y)i = Syi — y'lSxi — tj 



small as we choose by taking a sufficiently small value of Sxi. The 
substitution of this value for (0y)i in (2) gives 



S *'(-)-»' + - 



) +'!>(':, y, y') tei 
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and as (3) is to hold independently of the value of dxi, 

[|;»'M-»')+*fe!/,!/')]-0. (4) 

Comparing (4) with (2) we see that it is (2) with 5y replaced by 
d>l/(x) — dy, and Sxj by dx. 
If Xn is varied as well as x we have in addition to (2) the condi- 

and if y = fix) is to start at a point on the curve y = x(^) wo must 
replace 3y and Sxa by dxix) — dy and dx respectively. 

19. (a) Suppose that in our Minimum Surface problem of Art. 
7 the required generating curve is to be terminated by the point 
(0, j/o) and the given curve y = \('{x). 

We get as in Art, 7 

as the differential equation of the curve, and the condition 

\-^L=5y + y V 1 + y'= Sx,] = 0, 
LV 1 + y'^ i^^^ 

which by Art, 18 gives us 

|"-I|L= mx) -y')+v VT+7^1 = 
LVl + y'^ i = ^ 

or fy'^'^^^j + 1 = **■ 

Therefore the required curve meets the given curve at right angles. 

It has been shown in Art. 7 that it is a catenary having the axis 
of revolution as its directrix. 

(6) As a second example let it be required to find the curve of 
given length and terminated by the axes, which with the axes en- 
closes the maximum area. Here 



J=J[y + \^l + y'^dx. 
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57 -f\Sy + X ^^^^_^6y'jrf3; + |^y + X V 1 + j/'^j Sx, 



1 - 7 -Ji — . n 

'^^ V 1 + y'2 


(1) 


[ '"' 1 


(2) 


[^T%i<»~"''+'' + "^ +"'-]-»• 


(3) 


(1) integrates into 




ix + C)^ + (y - i:)' = X^ 





a circle. 

(2) gives [y'] = 0, and the circle crosses the axis of Y at right 
angles. 

(3) gives ■ = since y ~ when x = Xi, and the circle 

LVl + y'^^^ 

crossestheaxisofX at right angles. Therefore the area in question 
is the quadrant of a circle. 

20. The Principle of Least Action. 

If r is the kinetic energy of a moving system the time-integral of 
2T is called the action of the system. If 7 is the potential energy of 
the system and the forces are conservative it is well known that the 
so-called equation of energy T + Y = h, where A is a constant, holds 
good during the motion. 

The Principle of Least Action is that in the actual motion of the 
system the action is less than in any hypothetical motion, between 
the same initial and final configurations, in which the equation of 
energy holds good. 

In Hamilton's Principle the actual motion is compaj-ed with a 
hypothetical motion in which initial and final configurations and. 
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time of transit are the same as in the actual motion, and to bring 
about the hypothetical motion additional forces would have to be 
imposed on the system, and these forces would do work. 

In the Principle of Least Action the actual motion is compared 
with a hypothetical motion in which initial and final configurations 
are the same as in the actual motion, but as the equation of enei^y 
holds good, the time of transit is different in the two motions ; and 
the hypothetical motion could be brought about by imposing con- 
straints that would do no work. 

The formulation of the Principle of Least Action is 

S \2Tdi = 0, T -^V = h. 



As in the case of Hamilton's Principle we shall prove the Prin- 
ciple of Least Action for a single moving particle, and we shall use 
the notation and as far as possible the results of Art, 12. 

Keeping in mind that the potential energy is the negative of the 
force function, i.e. V — — TJ, and that d, the time of transit, is 
varied. 

bj2Tdt = sJ\t + h- V)dt = f(ST - SV)dt + [T + h- V]Sti. 
fsiT - V)dt = m[x'Sx + y'Sy + z'Ss], '' 



P^ 



as in Art. 12. 

Since the initial positions are the same in the two motions.. 
5x = Sy = 52 = when t = 0. Since the final positions are the 
same but are reached in the time d in the actual motion and in the 
time d -I- Sd in the hypothetical motion 5x = 5y — 5z = when 
t = ii+ Sti, and are approximately — x' 5d, — y' 5d, and — s' 6d 
when i = d- Therefore (v. Art. 18) 

f S(T - V)dt = - m[x'^ + y'" + 2'=]5d = - [2r]6d 
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approximately, and 

5j2Tdt = ~-[T+V - h]Sf, approximately. 
=0. 

Ih any concrete problem of course the Principle of Least Action 
leads to the same differential equations of motion as Hamilton's 
Principle, 

As an example we shall take the problem of planctarj' motion 
(jj.Art. 12). 

r - f (r'^ + r^'^), V ^ - ^. as in Alt. 12. 

By the Principle of Least Action 

SI = 5 jm(r'^ + rV'^)rf( = 



and 1(r" + r-i,"}-'^.h. 

"'"^/[^ ('" + ''*") +l> + l~^]dt 
= fm\r'Sr' + r'^'S^,' + r*'"Jr - ^^Srldl 

Integrating by parts 



+ m[r'5r + rV'S^t + [f (r'^ + r'<k") +h + ^1 St,. 
The terms outside the integral sign reduce to 
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and are equal to zero. Therefore 



as in Art. 12. 

And we get r" — r(f>'^ = — 

and rV' = ^, as in Art. 12. 

As Hamilton's Principle is more general than the older Principle 
. of Least Action and is decidedly easier to use it is much more com- 
monly employed in practice than the latter principle. 

21, The whole subject of Variation of the Limits can be treated 
from a different point of view, that of parametric representa- 
tion. 

Instead of tiying to determine y as a function of x so that 

/ = j 4>ix,y, y')dx shall be a maximum or a minunum let us 

start to determine x and y as functions of a parameter t. 

We shall find it convenient to represent — by x and ~ by y,- in 

which case y' = -, and dx — xdr. t of course will be our independent 

variable, and we shall suppose that m and n are fixed while the 
corresponding values xa and Xi may be varied. 

/ = / 4,{x, y, y') dx = { <t>{x, y, y')xdT. 
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Integrating by parts 



-/{l 


dx dry' *aWJ M % 






^-[(*-ii--l 


«»]• 


To make &I ■■ 


= we must make 






dy dr&y' 


(1) 




^^-^(^^ii- 


(2) 




(*-il>-+>.- 


(3) 


when r = m, 






and 


(*^ii-.+l*-« 


(4) 


when T = T,. 







Dividing through by x (1) and (2) become 



(5) 



8y dx %' 
(3) and (4) reduce to 



(8) 

when X = Xi. 

(5) is Lagrange's Equation, (6) wliich can be shown to have tho 
same solution as (5) we can disregard. (7) and (8), if we have a 
relation between Sxa and Sya, and between 5xi and Sj/i, serve to de- 
termine the two arbitrary constants in the solution of (5). 



y Google 



PRINCIPLE OF LEAST ACTION 45 

For instance if the required curve is to be drawn from a given 
curve y — x(^) to a given curve y — 'p{x) as in Art. 18, 
^Va = x'(^) 5^0 and dyi = yp'(x) Sxi, approximately, aa in Art. 18; 
and (7) and (8) reduce to 

+ (x'(x) - V') ^' = ^ when X = Xo 

and 4< + {'P'ix) — y') ■»— / = when x ^ Xi v. Art. 18. 

22. To illustrate the practical working of this method we shall 
apply it to the second example in Art. 19, to the proof of the Prin- 
ciple of Least Action (w. Art. 20) and to the problem solved in Art. 
20. 

(a) iv. Art. 19). 1 = f[y + X Vl + y'^^jdx. 

If X and y are functions of a parameter t,Sx = and 5y is arbitrary 
at the start, i. e. when t = m; and dx is arbitrary and Sy = at the 
finish, i. e. when r = ri. 



/=/[!, + > 


V^+lh 


--jfN + x> 


/i' + 


WVr- 


ai./[xi, + i,te + 'i<^»^-.+-M]* 






<'[*- 


dr ^Ja? + yA 


^-[* + .4; 




M'' 


-'\Jh^'H'-. 


/.r> + W J_ 






Shice a/ = 


d \y 


. 








dr V a? + y= 




y-\- 


d \x 


= 0: 
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J CALCULUS OF VARIATIONS 


■Ji' + y' 


wllRn r = To 


, >>* n 


when r = Ti. 


Vi' + j' 


These equations reduce to 




1 ~} ^»' -0, 




S' + ,^ ' 0; 




"' n 


when a^ = 


V 1 + j'- 


' -0 


when X = x,. 



(1) 

(2) 



(1) ia (1) Art. 18 and integrates into (x + Cy + (y ~ K)" = \\ 
as does (2). The curve is a circle; and as y' = when a; = and 
J/' = CO when y = 0, the circle eroaaes both axes at right angles and 
the centre must be at the origin. 

(b) (v. Art. 20). 

We shall suppose that the time, (, and therefore the coordinates x 
and y, depends upon an independent variable t in such a way that 
when r — Ta,t = Oin. the actual and in the hypothetical motion, 
and when r = n, i = d in the actual motion and t = t, + Sd in the 
hypothetical motion. Note that Sx and 8y are zero when t — to 
and when r = n, and that S( = when r = tc. 

C^Tdt - f{T + h- V)dt = f[T + h~ Vjidr. 

s C2Tdt = { [(ar ~ 5V)i + (r + a - v)bt]dT 



= JiiST ~ mi + 2T6(]dr. 
- hY = m{x"hx + y"hy + z'lz). 
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PRINCIPLE OF LEAST 


ACTION 






x"5x 


1 d/i,\ igi 










- SV 


m d xSx + yby + 
~ t'dr I 


iSz 


- l^m 






m d ■x5x + ySy + 
~ t dr i 


z6z 


-ST- 


-2T 


Si 
i 




^^^ d xbx-\-iihj + . 


i5z_ 









[ST - dV]i + 2TSi - 

5 1 2T5t = mlx'Sx + y'Sy + z'Bs]='o 

since Sx = &y = &z = when r = tq and when r = ri. 
(c) Planetary motion (ii. Art. 20). 
We are to malce 

I.f[lir" + r^'')+k + ^-^]d, 

a minimum, while keeping 

j(r= + r* ) --j^ -/I. 
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